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ABSTRACT 

We present results from two new near-infrared imaging surveys. One survey covers 
47.2 arcmin^ to K{3a) = 20"* whilst a second, deeper catalogue covers a sub-area of 
1.8 arcmin^ to K{3(t) = 22.75™. Over the entire area we have extremely deep optical 
photometry in four bandpasses {UBRI), allowing us to track the colour evolution of 
galaxies to very faint magnitude limits. 

Our K— band number counts are consistent with the predictions of non-evolving 
models with < qq < 0.5. The K— selected {B — K) galaxy colour distributions from 
our surveys move sharply bluewards fainter than K 20. At brighter magnitudes 
{K < 20™) our K— selected {B — K) distributions indicate a deficiency of red, early- 
type galaxies at z ~ 1 compared to the predictions of passively evolving models, which 
implies either a pure luminosity evolution (PLE) model where star-formation continues 
at a low level after an initial burst, or dynamical merging. At fainter magnitudes, the 
continuing bluewards trend observed in {B — K) can be explained purely in terms of 
passively evolving PLE models. We detect 0.5±0.1 galaxies arcmin"^ with (I — K) > 4 
and 19 < < 20 mag. Althoug h this is a factor of ~ 3 (2ct) more objects than 



the recent survey of Barger et al. this is still lower than the predictions of standard 
passively evolving models and more consistent with PLE models containing small 
amounts of on-going star-formation. Our observed numbers of (/ — K) > 4 galaxies 
at X ~ 20 are lower than the predictions of passively evolving models or PLE models 
with a low level of continuing star- formation, suggesting that at least part of the larger 
deficiency observed in {B — K) at fiT ~ 20 may be due to star-formation rather than 
dynamical merging. 

As we and others have noted, the number redshift distribution at 18 < if < 19 of 
recent, deep K— selected redshift surveys is well fitted by non-evolving models, and 
passively evolving models with a Salpeter or Scalo initial mass functions predict too 
many galaxies with z > 1. Dynamical merging is one possible solution to reduce the 
numbers of these galaxies but (as we have suggested previously) a dwarf-dominated 
IMF for early-type galaxies could offer an alternative explanation; we show here that 
such a model reproduces well the optical-infrared colour distributions and K— band 
galaxy counts. 

Key words: Galaxies : evolution -Galaxies: photometry. 



1 INTRODUCTION 



In ourj )revious papers (iMetcalfe et al. 199l| 



1995; 



Metcalfe et al 



Metcalfe et al. 1996| ) we have investigated the evolu- 



tion of the counts and colours of faint field galaxies, culmi- 



nating in extremely deep optical imaging of a repre sentative 
region of sky we call the "Herschel Deep Field" (Metcalfe 



et al. 1998). In this paper we present new near-infrared ob- 



servations of this field, and combine this data with our ex- 
isting ultra-deep optical imagery. Near-infrared observations 
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Table 1. Photometric Limits of the Herschel Deep Field. 

Filter U B R I K 

Limit (3(t) 26.8 27.9/28.2 26.3 25.6 20.0/22.75 

Area (arcmin^) 46.4 46.8/2.8 48.5 52.2 47.2/1.78 



have many advantages which are oft-quoted in the literature, 
among them, the insensitivity of the k— correction to type or 
star-formation history at intermediate redshifts is perhaps 
the most relevant in a cosmological setting. 

Previously we have considered our observations in 
the context of pure luminosity evolution (PLE) models. 
Recently, additional support for these models ha s come 



from f ully complete, b lue-selected, redshift surveys ( Cowie 
Songail i, fc Hu 1996 ) which have revealed the presence of 



an extended tail in the redshift distribution, partially resolv- 
ing the paradox that non-evolving galaxy formation mod- 
els could be used to fit the results o f the early faint red- 
shift su rveys ([Glazebrook et al. 1995h ; Broadhurst, Ellis, & 
Glazebiook 1992). Such models, as we have shown, are able 
to reproduce all the observable quantities of the faint field 
galaxy population (counts, colours, redshift distributions) at 
least for low fio Universes and within current observational 
uncertainties. High fio Universes can be accommodated by 
the model if we add an extra population of low luminos- 
ity galaxies with constant star-formation rates wh ich boost 
the counts at faint (B > 25™) magnitude levels (Metcalfe 
et al. IS 196; CJampos 1997). In the PLE model, galaxies form 



monolithically at high redshift. Changes in luminosity and 
colour after the formation event are modelled with stellar 
population synthesis models. K— selected galaxy samples, 
at least until K ~ 20, are dominated by early-type galax- 
ies; therefore, changes in optical-infrared galaxy colours are 
particularly sensitive to evolution of this galaxy class. By 
investigating near-infrared number counts and colour distri- 
butions we can potentially address questions concerning the 
formation and evolution of elliptical galaxies. 

The two surveys presented here are complementary. 
The "wide" survey covers a large, ~ 50 arcmin'^ area to 
a shallow, K{Za) = 20™, depth, whilst the "deep" survey 
covers a much smaller area of ~ 1.8 arcmin'^ to a much 
fainter 7^(30-) = 22.75™ limiting magnitude. In combina- 
tion, both surveys provide a K— limited sample extending 
from K to K = 22™. 75. 

This paper is organised as follows: Section ^ describes 
the data reduction techniques used in both the wide and the 
deep surveys; in Section g number counts and colour distri- 
butions for our both data sets are presented; in Section ^ we 
discuss these results in terms of our modelling procedure; 
and finally, in Section ^ we summarise the main conclusions 
we draw from our work. 



2 OBSERVATIONS AND DATA REDUCTIONS 



2.1 Roducing tho "Wide" Survey 



The observations were made over four nights in July 1995 
at the 3.8m United Kingdom Infrared Telescope (UKIRT) 
using the IRCAM3 near-infrared detector with a standard 
K— (2.2 ^m) filter. Conditions were photometric on two of 



the four nights. The mean seeing was ~ 1" FWHM. The 
centre of the region we imaged, RA Oh 19m 59s Dec -1-0° 
4' 0", corresponds to the centre of our de ep optical fields 
which are described in Metcalfe et al. (199^ ). For this region 
we have imaging data complete to the limiting magnitudes 
listed in Table fu The optical observations were made at 
the 4.2m William Herschel and the 2.5m Isaac Newton tele- 
scopes (WHT and INT hereafter) and cover approximately 
~ 50 arcmin^. 

The IRCAM3 detector, a 256 x 256 element InSb chip, 
has a pixel scale of 0.286" and consequently each image cov- 
ers only ~ 1.5 arcmin'^. For this reason we observed a mosaic 
of 6 X 6 frames which completely overlapped the previously- 
observed optical fields. The entire mosaic was observed in 
each pass for two minutes at each position, and offset slightly 
for each observation from the previous pointing. In total, 
the entire survey area was covered fifteen times, giving a 
total integration time at each sub-raster of 30 minutes. The 
short exposures for each image is a consequence of the high 
sky brightness in the near-infrared and is fact the longest 
exposure time possible with IRCAM3 in the K— before sky 
background causes the detector to saturate. Our data reduc- 
tion techni que for both infr ared surveys follows the methods 
outlined in pardner (1995). 



A two-pass procedure is adopted: for the first stage 
in the data reductions a dark frame is subtracted from 
each night's images, calculated from the median of the dark 
frames taken that particular night. Next, the images are 
sorted in order of time of exposure. Then, a running sky 
is constructed from the median of the six images nearest 
in time to the frame being processed. This is necessary as 
the sky background varies on exceedingly short time-scales 
(i.e., under 30 minutes). This running sky is then subtracted 
from the current image and the procedure repeated for all 
the images in the mosaic. Next, for each night a "superflat" 
is constructed by medianing together all the observations for 
that night. Each image is then divided by the superflat. 

Once all the images have been dark- and sky-subtracted 
and flat-fielded they are grouped into 36 separate pointings 
and each stack averaged together using the IRAF imcom- 
BINE task. Next, masks were constructed from each of the 
stacked images by flagging all pixels more than 2>a from 
the sky background. These masks were then used in re- 
calculating the running sky frames. This was necessary as 
bright sources present in each of the individual exposures 
caused raised counts in the original running sky frames. The 
whole procedure outlined above was then repeated, using 
these new frames. 

Finally, the 36 stacked sub-rasters were combined into 
a single image by matching each sub-raster with objects on 
the extremely deep band frame and using this to cal- 
culate the correct position of each sub-raster relative to the 
others. The final image, after trimming, covers an area of 
47.2 arcmin^, overlapping our WHT deep optical fields. 
We calibrated this data by repeate dly (~ 30 times) ob- 



serving the faint UKIRT Standards (Casali & Hawarden 



1992) each night at a range of airmasses, although fortu- 



itously the mean airmass of the standards and of the pro- 
gram field images was approximately the same, so accurate 
determination of the extinction coefficient was not neces- 
sary. The rms scatter between our standards measurements 
was ±0.03 magnitudes. 
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2.2 Reducing the "Deep" Survey 

Observations for the "Deep" survey were made over six 
nights in 1994 September at UKIRT using IRCAM3 with 



a K-^ baud fllLyi. CundiLiuns wyiy pliuLumyLiic duiing Lhe 
first tw o niglita, and tlic rcmaiudcr of the data were acalcd to 
the photometry of bright objects in the field made on those 



nights 



ts. l iUe seeing ranged hum 0.8^^ Lu lA" . The LuLdl ub- 
serving himc on the field waa juat over 100 kaco (27.9 houra), 
and was made up of 838 individual exposures of 120 seconds 
each. Each individual exposure was made up of 12 co-adds 
of 10 seconds each. The telescope was dithered using the 
offset guide camera in a diamond pattern of 13 positions, 
each position separated by 5.25" in both Right Ascension 
and decUnation. Data reducti on p roceeded generally in the 
manner described in Section Dark frames were made 
several times each night, and were subtracted from the im- 
ages. A running sky was constructed, using the median of 12 
images (six before and six after). The running sky was sub- 
tracted (after applying the two-pass object masking system 
described above), and a "superflat" was constructed for each 
night and divided into the data. The input dither pattern 
was recovered and adjusted by a centroid of the brightest 
object in the field. The data taken each night were shifted 
and co-added. The photometry of each night's data were 
scaled to that of the photometric first two nights (which 
were themselves consistent). The final image covers an area 
of 1.78 arcmin'^, and has a limiting depth of 3cr — 22.75. The 
deep field is centred at 00 19 59.6 -1-00 02 56 (1950.0) and 
covers 81" square (140x140) 0.58" pixels. 



2.3 Object Detection And Photometry 

For both the deep and wide images we use the same prescrip- 
tion for object detection and photometry as in our optical 



studies ( Metcalfe et al. 1995 



Metcalfe et al. 199S ). Note that 



to improve the reliability of image detection both images 
were binned 2x2 (the 0.286" pixel size of IRCAM3, com- 
bined with median seeing of ~ 1" FWHM, means that our 
raw images are oversampled) . We fit the sky background ap- 
proximately with a 3rd order polynomial and then subtract 
this from the frame. Deep isophotal image detection is then 
undertaken (to a magnitude a bou t 1™ deeper than the 3(t 
limit for the frame - see Table 2.2) and the objects removed 
from the frame and replaced by a local sky value. Q The re- 
sulting image is heavily smoothed (with several passes with 
box smoothing up to 10 pixels on a side) and then subtracted 
from the original. The result is an extremely flat image on 
which the isophotal detection is then rerun. At this stage, 
in order to reduce false detections, any images within a few 
pixels of on e another ar e recombined into one. A Kron-type 
magnitude ( Kron 198C| ) is then calculated for each image 
using a local sky value. This is essentially an aperture mag- 
nitude calculated to mRk, where Rt is the Kron radius and 



m is the multiplying factor given in Table 2.2. It is necessary 
to set a minimum aperture equivalent to that for an unre- 
solved object, and apply a correction to total magnitude. 
These are also listed in Table 2.2. A more detailed discus- 



sion of this technique is given in Metcalfe et al. (1991 ). Our 



* Our 3a- limit is defined as when the varianee of the sum of the 
pixel values within our detection aperture reaches this value. 



optical-infrared colours are measured inside fixed 1.5" radius 
apertures. Note that to form the "Deep" (B — K) colours, 
the K— band magnitudes are matched with the 2' diameter, 
stacked, ultra-deep B— band image described in Metcalfe 



et al. (1998) (and see Table pi). All other colours are formed 



from matches with the 7' x 7' Herschel Deep Field optical 



image s. Star-galaxy separation is taken from Metcalfe et al 
(1998D. 



2.4 Confusion corrections 

The deeper the ground-based number-counts are extended 
the more important the corrections for confusion becomes. 
This is of particular relevance for the deep K-band data. 
To make an estimate of the completeness of our counts on 
this frame we have added numerous artificial stars of various 
known magnitudes to the real data frame and then reanal- 
ysed the resulting images using the normal data reduction 
procedure. We then compare the measured magnitudes with 
those which we input. Table 2.4 gives the mean magnitudes, 
scatter, and detection rate for these stars. Note that an im- 
age is considered undetected if it is merged with another 
image and the combined brightness is a factor two or more 
greater than the true magnitude, or if it is not found within 
±2 pixels of its true position. As expected the detection rates 
in the real data drop as the magnitude falls. This is almost 
entirely due to objects being merged with other, brighter 
objects - isolated objects are recovered with almost 100% 
efficiency. Note, however, that the mean recovered magni- 
tudes are close to the true values. Although most objects at 
the limit of our data are galaxies, they are in general close 
to being unresolved, and so we correct our counts by the 
detection rates implied by our artificial stars. 

Due to the much wider mean separation between ob- 
jects, the shallower K— band data do not suffer from con- 
fusion, and no correction is made. 



3 RESULTS 

The final "wide" catalogue contains 298 galaxies to a 
K(io) = 20.0"* over an area of 47.2 arcmin^. The deep 
catalogue contains 86 galaxies to a Kifia) — 22.75™ over an 
area of 1.8 arcmin^. 



3.1 K— band Galaxy Number Counts 

Number counts derived from both our catalogues are pre- 
sented in Figure |l|, which shows our data (with error bars 
calculated using Poisson counting statistics) and a compi- 
lation of published number counts. At bright magnitudes 
(le"" < iC™ < 20™) we agree weU with the counts from the 
literature, whereas at fainter magnitudes (20™ < K < 22™) 



our data favours the higher counts found by Bershady, 



Lowenthal, fc Koo (1998 ) as opposed to the lower counts 



measured by Djorgovski et al. (1995). Bershady, Lowenthal 



fc Koo claim that the aperture corrections employed by 
Djorgovski et al. (1995) resulted in an overestimation of the 

our results support 



0.5" 



depths of the latter survey by 
this interpretation. 

It is worth not ing that in our last bin, as in Bershady, 
Lowenthal, fc Koc, our counts are only approximately 50% 
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Table 2. Detection and measurement parameters used in both surveys. 

Frame Limiting isopliote Detection limit Min. r Kron multiplying Correction to 

(mag/arcsec^) (K mag) (") factor total (mag) 

Deep K 24.50 23.50 1.25 1.40 0.32 

Wide K 21.50 22.00 1.35 1.50 0.25 



Table 3. Results of adding artificial stars to the deep K data. 

True magnitude Measured magnitude Detection rate (%) 



20.75 20.70 ±0.17 90 

21.25 21.22 ±0.21 88 

21.75 21.76 ±0.26 85 

22.25 22.36 ±0.36 76 

22.75 22.83 ±0.37 51 



comple te and therefore subiect to large corrections. Fur- ([Griffiths et al. 1994| and the Hubble Deep Field ( [Williams| 



thermo -e, given the small areas of surveys at this magni- 
tude level (where areas covered are typically ~ 1 arcmin ) 
it is als o possible that cosmic variance could account for the 



et al. 1996) indicate that counts for elliptical galaxies at 



large o bserved ficld-to-ficld variations between the different 
groups. 



intermediate magnitude levels (IS™ < /siew < 22™) fol- 
low t he predictions of a non-evolving model ([Driver et al 



1998 ; Driver, Windhorst, & Griffiths 1995 



Considering the slope of our number counts, we derive 
a value of d{\og N) / dm = 0.37 ± 0.03 (la errors) over the 
magnitude range covered by both wide and deep catalogues 
(namely, K ~ 16—23). This value is ir i good agreement with 



d{\og N) / dm = 0.36 ± 0.02 found by Bershady, Lowenthal^ 
do not find the count slope for the deeper survey 
to be significantly different for the shallower one (howe ver, 
as pointed out by Gardner, Cowie, & Wainscoat (1993) the 
compiled K— band number counts show a slope change at 
K ~ 17). 



In 



Figure [1[ we hIho plot the predictions of two non- 
evolvin g galaxy count models, one with deceleiation 
ter qo = 0.5 (solid line) and one with go = 0.05 (dotted line). 
Both models give a good fit the observations from isT ~ 14 to 
if ~ 22. These n on-evolving model were co mputed from the 
k— corrections of Bruzual & Chariot (1993) and the z = lu- 
minosity function parameters given in Table ^. We compute 
by combining the tabulated Mg with the rest-frame 
[B — K) colours. In our models we adopt the same normal- 
isation as in our previous papers. As we and ot hers have 
noted before (Metcalfe et al. 1991; ^hanks 1989 ), normal- 
ising galaxy counts at brighter {B ~ 16™) magnitudes re- 
quires that substantial amounts of B— band evolution must 
take place at relatively low {z < 0.1) redshift. Conversely, by 
normalising at 18™ < B < 20™ we allow a better fit at in- 
termediate magnitudes and reduce the amount of evolution 
required to fit the counts at fainter magnitudes, but under- 
predict counts at brighter magnitudes. Support for this high 
normalisation ha s recently come f rom the 2-micron all-sky 
survey (2MASS; Kleinmann 1992 ) which finds a significant 
underdensity in the galaxy counts at 10™ < K < 13™ and 
is discussed further in a forthcoming paper (Shanks et al 
1999). 

Could substantial galaxy luminosity evolution be oc- 
curring at these low redshifts? Bright, K < 20™ sam- 
ples are dominated by early-type galaxies Morphologically- 
segregated number counts from the Medium Deep Survey 



Glazebrook et al 



1995a[ ). This would argue against an evolutionary explana- 
tion for the steep slope observed in the range 10™ < K < 
14™. We are then faced with finding an alternate explana- 
tion for the number cou nts under- density; one suggestion 
com es for the findings of Metcalfe, Fong, fc Shanks (199f: ) 



and Berlin & Dennefeld (1997) who both show how scale 
errors (i.e., non-linearities) in the photographic APM pho- 
tometry of Maddox et al. (1990) could result in anoma- 
lously low galaxy counts at brighter magnitudes in B— se- 
lected surveys; however, the K— band surveys are carried 
out with electronic detectors, a nd such effects are likely to 



be les s significant (although the Mobasher, Ellis, fc Sharplei _ 
(1986| ) survey was based on objects selected from B~ band 



photographic plates). Another suggestion is that our galaxy 
may reside in a locally under-dense region of the Universe 
(^hanks 1989); however, to produce the observed discrep- 



ancy in the number counts would require that this void 
would be ~ 150/i~^ Mpc, and have an under-density of 
~ 30%; fiuctuations on this scale are difficult to understand 
in terms of current cosmological models and also measured 
large-scale bulk flows. Some evidence for a local void has 
come f rom studies of pe culiar velocities of type la super- 
novae (Zehavi et al. 1995) although the size of the void they 



detect (~ 70h~ Mpc) is not large enough to explain the 
low number counts. For time being, resolution of this issue 
awaits a better determination of the bright end of the B— 
and K— number counts by much larger surveys. 



3.2 Optical-Infrared Colours 

Figure ^ shows the {B — K) vs K colour- magnitude diagram 
for both samples for galaxies in both the "wide" (open tri- 
angles) and the "deep" survey (open circles). Additionally 
the median {B — K) colour computed in one- magnitude bins 
is shown (filled squares) and compared t o the values given 
in Gardner, Cowie, & Wainscoat (1993) (filled triangles). 
Error bars on the median colours were calculated using a 
bootstrap resampling technique. 
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Figure 1. Differential galaxy number counts in half-magnitude intervals as a function of K— limiting magnitude as measured from our 
UKIRT wide survey (filled circles) and from our UKIRT deep survey (filled squares). Also plotted is a compilation of published K— band 
number counts and the predictions of the two non-evolving models described in the text for go = 0.05 (dotted line) and go = 0.5 (solid 
line). Error bars are calculated from the number of galaxies in each magnitude bin using Poisson y/n counting statistics. The inset shows 
the number counts from our surveys alone (this time with open symbols) compared to the non-evolving predictions. 



Incompleteness is represented by the dashed line; ob- 
jects to the right of this line have B > 28™. Several objects 
have been identified in both surveys which are undetected in 
B— but are detected in K— and these are plotted in the fig- 



ure as right-pointing arrows and are discussed in Section 4.3 
Considering the median [B — K) colour we see that it in- 
creases steadily until around K ~ 18™ where it reaches a 
maximum value of ~ 6. After this, it sharply becomes bluer 
and this blueward trend continues to the faintest limits we 
have measured. Our median colours agree well with those 



in Glardner, Cowie, & Wainscoat (1993), even for brighter 



bins where the number of galaxies in our sample are much 
smaller. We have also found our K < li [B — K) distri- 



bu tions to be in good agreement with the brighter survey 
of Saracco et al. (1997). Finally, it is worth noting the ex- 
tremely large spread in colour in this diagram - at if ~ 20 
[B — K) ranges from 2 to 9. 

Additionally, the predictions of the two non-evolving 
models shown in Figure |l] are plotted. In general, the go = 
0.5 model gives a slightly better fit to the observed colours, 
at least until K ~ 20™; at fainter magnitudes, galaxies are 
significantly bluer than both model predictions. 

Another interesting feature in Figure ^ is the apparent 
deficit of galaxies with magnitudes in the range > 21™ 
and colours from {B — K) > 5. The different areal coverages 
and depths of the two surveys doubtless exacerbates this ef- 
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Figure 2. K— magnitude against {B — K) colour for galaxies in the wide survey (open triangles) and the deep survey (open circles). Also 
plotted are the median (B — K) colours in one- magnitude bins for our survey (filled squares), with la error bars calculated by a bootstrap 
resampling technique. Triangles show median colours from the compilation of Gardner et al. 1993. The dotted line represents the region 
of incompleteness; galaxies below this line have B < 28™. Drop-outs are plotted at the upper limit of their colours as right-pointing 
arrows. For comparison, predictions from two non-evolving models are plotted; one for go = 0.05 (dotted line) and one for go = 0.5 (solid 
line) . 



feet but its presenee is consistent with a population of galax- 
ies which turns rapidly blueward faintwards of K ~ 20. We 
discuss the origin of this feature in Section 4.2. 

In Figure I we plot {B - R) vs. {R - K) colour for 
all the galaxies in the wide survey (open circles) and the 
deep survey (filled circles) . Fainter galaxies are plotted with 
smaller symbols. The "deep", K < 22.75™ sample is much 
bluer in both {B — R) and {R — K) than the "wide" than 
the brighter, K < 20™, sample, which reaches {R — K) ~ 6. 
We also plot colour-colour tracks representing the difi'erent 
galaxy classes in our evolutionary models. We defer discus- 
sion of these tracks until the following section. 



4 INTERPRETATION AND MODELLING 
4.1 Outline of the Models 

To interpret our results we investigate variants of a pure 
luminosity evolution (PLE) model in which star-formation 
increases exponentially with look-back time. Earlier ver sions 



In this paper we consider cosmologies in which A = and 
take Ho — 50 kms~^Mpc~^, although changing the value 
of Hq does not affect any of the conclusions of this paper. 
Two values of the deceleration parameter, go, go = 0.05 
and go = 0.5, are adopted, corresponding to open and fiat 
cosmologies respectively. (Zero-curvature cosmologies with 
Q. = 0.3 and Ao = 0.7 are almost identical to our go = 0.05 
model because co-moving distance and look-back time as 
a function of redshift are almost identical t o that found in 
the zero-Lambda models, at least to z ~ 2 ( Peebles J_984 ) ) . 



of these models are discussed in our previous papers (Met- 
calfe et al. iggi] ; [Metcalfe et al. 1995| ; [Metcalfe et al. 1996D 



The input parameters to our models (given in Table [^ con- 
sist of observed local galaxy parameters (namely, rest-frame 
colours and luminosity functions) for each of the five mor- 
phological types (E/SO, Sab, Sbc, Scd and Sdm) we consider 
in our models. These morphological types are divided into 
elliptical (E/SO) and spiral (the remainder) and these two 
classes are each given a separate star-formation history. We 
could, in principle, sub-divide the spirals into different mor- 
phological types each with different star formation histories 
but for simplicity we do not; (k + e) corrections for the differ- 
ent types are fairly similar to each other in these models in 
any case . Instead, taking a Sbc model as representative of 
all types we produce the other types by normalising the Sbc 
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Figure 3. (B — R) colour against (H — K) colour for galaxies in the wide survey (open circles) and the deep survey (filled circles). Fainter 
galaxies are plotted with smaller symbols. Also shown are evolutionary tracks for the five classes of galaxy in our model, as computed 
from the population synthesis code of Bruzual & Chariot (93). The solid line represents type E/SO, the dotted line Sab and the dashed, 
long dashed and dash-dotted lines types Sbc, Scd and Sdm respectively. The dot-dashed line shows the path of of an elliptical galaxy 
which has an IMF slope, x, of 3. For the elliptical tracks, the labelled stars indicate the redshift. 



track to the observed rest-frame (B — K) co lour s from Table 



^ As we have already explained in Section 3JJ our (/>' s are 
chosen to match the galaxy counts at B ~ 18 — 20 and we 
seelc to explain the low number counts at bright magnitudes 
from a combination of photometric errors and anomalous 
galaxy clustering, rather than substantial evolution at low 
redshift. Our models also include the effects of the Lyman-a 
forest, and, for spiral types, dust extinction c orrespondi ng 
to the Large Magellanic Cloud as described in Pei (1992). 



PLE models have difficulty in correctly reproducing the 
observed K— band redshift distribution. Even essentially 
passive evolution over-predicts the numbers of galaxies seen 
at 1 < 2 < 2, and this fact has been used by s everal au- 
thors (most recently Kauffmann & Chariot (1998)) to argue 



for the hierarchical, merger-driven galaxy formation mod- 
els in order to reduce the mean redshift of K— selected 



galaxy redshift distributions. However, in Metcalfe et aL^ 
(1996n~]we showed how PLE models could produce a K— 
ba nd redshift distribution compatible with the observations 



of |Cow]e, Songaila, fc Hu (1996|) by assuming a very steep 
slope {x = 3) for the stellar initial mass function (IMF) in 
E/SO types. This steep, low-mass-star dominated IMF re- 
duces the amount of passive evolution observed in K— and 



allows us to reproduce the observed redshift distribution, 
which is close to the predictions of the non-evolving model. 
It s hould be noted however, tha t if the incompleteness in 
the [Cowie, Songaila, fc Hu (1996 ) n(z) were entirely due to 
z > \ galaxies then the deficit relative to the PLE model 
(with standard IMF) is only ~ Scr from \fN statistics alone 
and is likely to be less significant when galaxy clustering is 
taken into account. Conversely, the B— band redshift dis- 
tribution from the same survey does have a large fraction 
of high redshift (z> 1), high luminosity {L > L") galaxies. 
However, this is not a problem for our models. Our steep 
luminosity function slope for spiral types combined with a 
small amount of internal dust extinction (corresponding to 
Ab ~ 0.3™ at 2 = 0) allows us to reproduce this result. 

Although the resolution of the dataset discussed in this 
paper is too low to extract useful size information for the 
faint galaxy population, we mention in passing that observed 
faint galaxy sizes are compatible with these m odels, models. 



as w e have demonstrated in a recent paper (shanks et al 



199^ ). We have used a modified version of Howard Yee's 
PPP code to produce a simulated HDF frame which con- 
tains galaxies with parameters based on the output of our 
PLE models. For this image, spiral galaxies were generated 




Figure 4. Predicted redshift distributions for botli "Scalo" (t = 
1.0 Gyr) and "x=3" (r = 2.5 Gyr) models (dotted lines and solid 
lines) compared with the results of the K— (lower panel) and B— 
(upper panel) redshift distributions taken from Table 1 of Cowie 
et al. 1996. Incompleteness is represented by the open boxes. Also 
plotted are the predicted redshift distributions for elliptical types 
(long dashed and short dashed lines) for the "Scalo" and "x=3" 
models respectively. 



assuming that luminosity evolution only affects disk s urfac e 
brightness and not size, combined with Freeman's (1970), 
law to relate disk size to absolute magnitude at z=0. For 
bulges, we adopted the the d iameter-magnitude relation of 
Sandage & Perelmuter (199C). By applying our photometry 



package to this image we account for selection effects in a 
realistic manner. We find that the galaxies in the simulation 
have sizes not dissimilar from those measured in the HDF, 
in part because many are intrinsically small galaxies seen 
well down the LF. 

In Figure ^ the effect which changing the IMF slope 
for early types has on predicted redshift distributions is il- 
lustrated. The dotted line shows a mo del computed using 
T = l.Q Gyr and a Scalo ([Scalo 1986|) IMF, whereas the 



solid line shows the x = 3 model which has r = 2.5 Gyr. Also 
plotted are the redshift distri butions taken from Table 1 of 
Cowie, Songaila, fc Hu (199^ ). Objects which were marked 



model predictions to the total number of observed objects 
(i.e., unidentified objects are included). We take go = 0.05 
and Zf — 6.4 in all our models. The star-formation e-folding 
time T, is 9.0 Gyr for all the spiral types. The a:: = 3 model, as 
it dominated by dwarf stars produces much smaller amounts 
of passive evolution in the K— band and as a consequence 
the median redshift of galaxies in this model is much lower 
than in the Scalo model, which produces an extended tail in 
the redshift distribution. Notice the effect of changing the 
IMF slope in the B~ band is negligible. 

This point is further illustrated in Figure ^ where we 
show the median redshift of elliptical galaxies with a range of 
star-formation histories for several K— selected magnitude 
slices. This panel shows that, in all all magnitude ranges, el- 
liptical galaxies in the "Scalo" model have a higher median 
redshift than in the "x=3" model, irrespective of r (demon- 
strating that the differences in ^ are due to changing the IMF 
and not r). This point is discussed further in Section 4.2. 

In Figure ^ we illustrate the effect which varying the 
IMF slope and r has on the E/SO colour-redshift relation. 
The solid, dotted and long-dashed lines shows Scalo IMF 
tracks with r = 0.5 Gyr, r = 1.0 Gyr and r — 2.5 Gyr 
respectively, whereas the long dashed lines shows an a; = 3 
model with r — 2.5 Gyr. As is apparent from the plot, at 
z ~ 1 {B — K) colour depends very sensitively on the as- 
sumed value of T and the slope of the IMF. For a given IMF, 
longer r's cause the peak in the {B — K)-z relation to shift 
to progress ivel y higher redshifts, a point we will return to 
Finally, the x 



in Section 4.2 



as unidentified in each band in this table are shown in their 
respective incompleteness boxes. We have normalised both 



3 IMF model track, as it 
is dominated by low-mass stars, is in general redder (for a 
given value of r) than the Scalo IMF track. We note that the 
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Table 4. Luminosity function parameters and rest-frame colours used in all models 



Parameter / Type 


E/SO 


Sab 


Sbc 


Scd 


Sdm 


i>* {Mpc-3 X 10-*) 


10.2 


5.09 


6.82 


3.0 


1.5 


a 


-0.7 


-0.7 


-1.1 


-1.1 


-1.5 


M% 


-21.0 


-21.0 


-21.32 


-21.44 


-21.45 


(B^K) 


3.93 


3.78 


3.51 


2.90 


2.26 


{I-K) 


1.57 


1.63 


1.7 


1.47 


1.0 



^ 4 
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Figure 6. a: (lower panel) (B — K) colour against redshift, z, 
for an elliptical galaxy with an x = 3 IMF and r = 2.5 Gyr 
(solid line), and for the Scalo IMF and with three values of r; 
T = 0.5 Gyr (dot-dashed line), r = 1.0 Gyr (dotted line) and 
T = 2.5 Gyr (long dashed line). Also shown is a non-evolving 
model (short dashed line). b:(upper panel) As for the lower panel, 
but with {I — K) colour. 



Figure 7. Galaxy number counts for two evolutionary models - 
the "scalo" model, (solid and dotted lines), and the "x=3" model, 
(short dashed and long-dashed), for low and high values of qg 
respectively. K— band counts are shown in the upper panel; B— 
band ones in the lower one. 



single-burst models of Worthey (1994) have shown that the 
amount of passive evolution in the rest frame /— and R— 
bands is sensitive to metallicity variations; consequently at 
z = 1 — 2 the predicted K— magnitude from our essentially 
single-burst model is also expected to be approximately in- 
dependent of metallicity. Therefore reasonable changes in 
metallicity seem unlikely to change the rate of passive evo- 
lution at K— which, in the context of these models, only 
appears to depend on the slope of the IMF." 



We now return to Figure q where we have over-plotted 
on our (B — i?) — (i? — K) diagram evolutionary tracks for 
each of the galaxy types in our model. The dot-dashed and 
solid lines shows the path of an elliptical, whereas the re- 
maining tracks show the four spiral types in the models. On 
the elliptical tracks, the labelled stars indicate the redshift 
at selected intervals. Generally, there is good agreement be- 
tween the model tracks and the observed colours. The low 
redshift (0.5 < z < 0.0) elliptical track is populated by 
the brightest galaxies in our sample. Generally, we do not 
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find a large fraction of objects have colours which cannot be 
reproduced by the trac ks (in contrast with the findings of 
Moustakas et al. (1997) who claimed there was a significant 
population of objects whose colours could not be reproduced 
by the models). There are some objects which appear to be 
redder than the tracks, but these do not constitute a signif- 
icant fraction of the observed galaxy population. The broad 
range of {B — K) colours of objects with {R~ K) > 4 (which 
according to the models should all be E/SO galaxies) is a real 
effect which we interpret as evidence of the diverse range of 
star-formation and metallicities histories present in the el- 
liptical population. 

4.2 Counts and Colour Distributions 

In Figure Q we present a compilation of all published galaxy 
number counts in B— and K— bandpasses, as well as the 
predictions from our two evolutionary models. 

What is immediately apparent is that both the mod- 
els discussed provide a good fit to the number counts over 
a wide magnitude range for the low qo case. Furthermore, 
in both B— and K— bandpasses at faint magnitudes the 
differences between the Scalo and "x=3" models is much 
smaller than the differences between the high- and low-go 
mod els. It is also apparent, as has been noted by many oth - 



Babul fc Rees 1992) 



ers ( [Metcalfe et al. 1996^ [Campos 1997 
that, unless another population (such as bursting dwarfs) is 
added to the models then the qo ~ 0.5 model under-predicts 
the observed B— band galaxy counts. In the K— band this 
model is marginally excluded at the faintest magnitudes. 

Number counts are a coarse test of any model; the ex- 
tremely deep optical photometry covering the Herschel Field 
is one of main virtues of this work. So, to test our models 
more stringently we now investigate our observed optical- 
infrared colour distributions. We confine our discussion to 
the qo = 0.05 case as these models provide a better fit to 
the counts, although the colour distribution for the go = 0.5 
models will be very similar. 

In Figure ^ we compare in detail our model predictions 
with our observations for K— selected {B — K) colour dis- 
tributions and compare these with the data (shown as the 
solid histograms) in six one-magnitude slices from K — 16™ 
to K — 22™. The "Scalo" model is shown as a dotted and 
a dot-dashed line and the "x=3" model is represented by a 
solid line. All models have been convolved with a conserva- 
tive 0.3™ colour measurement error (overestimating errors in 
the brighter magnitude slices). For the non-evolving model, 
as for the evolving models, fc— corrections are calculated 
from the models of 



Bruzual fc Chariot (1993) 



Incompleteness in the histograms was determined by 
counting how many galaxies had K— magnitudes but were 
not detected in B— in each magnitude slice (and conse- 
quently have undetermined colours). The colour of the red- 
dest object which could be detected is represented by the 
downwards-pointing arrow, which is plotted at the colour 
corresponding to the K— magnitude at the centre of each 
magnitude slice (these objects have B > 28™). The num- 
ber of non-detections is shown by the size of the box. Mod- 
els have been normalised to the total number of galaxies in 
each slice; in the latter panels this includes B— band non- 
detections. 

The non-evolving model predictions (dashed line) 



rapidly diverge from the data. Faintwards of the K = 
17™ — 18™ magnitude slice a bump appears in the model 
colour distributions, corresponding to the unevolved, red, 
{B — K) ~ 8 elliptical population. At approximately the 
same magnitude limit and at bluer colours (this time at 
{B — K) ~ 4) a second peak becomes apparent. This corre- 
sponds to the model spiral population, and at fainter mag- 
nitudes this peak becomes the more prominent of the two. 
This is a consequence of our adopted luminosity function 
parameters and morphological mix. 

Qualitatively, the evolutionary models follow the pat- 
tern of the non-evolving model, though the location and am- 
plitude of the E/SO peaks depend sensitively on the choice 
of IMF and r. On the blueward side of the distributions, 
the spiral model colours appear to be a good fit to the ob- 
servations. Referring to Figure ^ it is apparent that beyond 
z ~ 0.3, the {B - K) colour for the r = 0.5 Gyr E/SO 
model is redder than the r = 1.0 Gyr model; furthermore, 
as the median redshift of both populations is almost the 
same (Figure |5|) , changes in the colour distributions are al- 
most completely a p roduct of differences in th e {B — K)-z 
relation. In general, Bruzual fc Chariot (1993 )-type models 
for elliptical galaxies with short r's predict a (B — K) — z 
relation in which galaxies are at their reddest at z; ~ 1, and 
become bluer thereafter until z ~ 3 at which point they turn 
redder once more as a result of Lyman-a forest reddening. 
The position of this turn-over depends on the assumed r, 
the IMF and the redshift of formation zf . 

If we now we consider the "x=3", r = 2.5 Gyr model, 
we see that its colour distributions are in general bluer than 
those predicted from the "Scalo" r = 0.5 or r = 1.0 Gyr 
models. This is partly a consequence of the longer r adopted 
in this model, but also of the fact that (as we have seen in 
Figure ^) galaxies lie at a lower redshift than the "Scalo" 
model and consequently are sampling a bluer region of the 
(B — K)-z relation. So, despite the x — 3 model being low- 
mass star dominated, the larger value of r allows us to pro- 
duce colours which are blue enough to match the observa- 
tions. 

Returning to Figure ^ we can now understand the ori- 
gin of the differences between the {B — K) colour distribu- 
tions predicted by the models. Considering first the brighter 
magnitude slices in Figure |, at /sT = 17.0™ - 18.0™ we 
see that the "Scalo" r = 1.0 Gyr model predicts a peak 
at {B — K) ~ 7 which appears in the "x=3" model at 
{B — K) ~ 6. At this magnitude, elliptical galaxies in the 
"Scalo" model have Zmcd ~ 0.6, whereas for the "x=3" 
model Zyn^d ~ 0.4 (from Figure^; as a consequence "Scalo" 
model colours are redder than "x=3" ones. Furthermore, the 
essentially passively-evolving Scalo r = 0.5 Gyr model con- 
tains insufficient recent star-formation to shift the peak at 
[B — K) ~ 9 blueward; as we have seen from our considera- 
tions of the non-evolving model, this peak is from unevolved 
ellipticals. At fainter {K > 20™) magnitudes, both the x = 3 
model and the Scalo r = 1.0 Gyr model turn sharply blue- 
ward; the Scalo r — 0.5 Gyr model colour distribution, how- 
ever, has a extended tail which reaches {B — K) ~ 9 — 10. 

It is clear from Figure ^ that either a Scalo IMF with 
a, T = 1.0 Gyr or an a; = 3 IMF with t = 2.5 Gyr model re- 
produces the main features of the data histograms, including 
the broadening of the colour distribution at the K < 18™ 
slice, and its subsequent narrowing and blueward trend in 
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Figure S. K— selected {B — K) colour distributions (solid histogram) in six magnitude slices for both the wide (brightest four slices) and 
the deep survey (faintest two slices). Also shown are predictions from four models: "x=3" (solid line), "Scalo" (dotted and dot-dashed 
lines) and "non-evolving" (long dashed line). The two "Scalo" models have differing star-formation histories; the dotted and heavy-dashed 
line has r = 1.0 Gyr whilst the heavy dot-dashed line has r = 0.5 Gyr. The x = 3 models has r = 2.5 Gyr. The model predictions have 
been normalised to the total number of galaxies in each slice. In the fainter {K > 18™) bins this includes the B— band non-detections 
discussed in Section 4.3 and the number of such objects in each magnitude slice is represented by the incompleteness boxes. The position 
of the downward-pointing arrow represents the colour of the reddest object at the centre of each magnitude slice which could be detected; 
objects in the incompleteness boxes must lie rightwards of this arrow. For all models in this plot we adopt go = 0.05. 



the fainter slices. Furthermore, as we will see in the following 
Section, these models are consistent with the redwards lim- 
its of our data. Of these two mod els, t he x = 2>, t = 2.5 Gyr 
model as commented in Section 4.1, gives a much better 
fit to the observed K— band redshift distributions. As is 
apparent from Figures I and I distinguishing between the 
effects of changes in IMF slope or the star-formation his- 
tory on the colour-redshift relation is difficult. Making this 
choice requires additional information such as the redshift 
distribution of K— selected galaxies. 



Finally, these models also allow us to understand the 
deficit of galaxies in the colour-magnitude relation at {B — 
K) > 5 and AT > 21™ which we commented upon in Sec- 



tion 



3.2. Firstly, as a consequence of our choice of luminosity 
function parameters elliptical counts turn over at if ~ 20 
and spiral types dominate faintwards of this. Secondly, from 
Figure ^ we see that at if ~ 21 E/SO's are at z ~ 1; from 
these redshifts to 2 > 3 Figure m indicates that elliptical 
colours turn sharply blueward. In combination, these fac- 
tors produce the rapid blueward trend observed at if > 21™ 
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Figure 9. K— selected (/ — K) colour distributions for both the wide and deep surveys (solid histogram). All symbols are as in Figure [ 



and the absence of redder objects, which happen in all our 
models regardless of our choice of IMF or r. 

4.3 (/ — K) distributions and extremely red 
galaxies 

In the previous section we discussed the broad character- 
istics of the K— selected {B — K) colour distribution and 
showed how these features can be understood in terms of 
our models; in this section we turn our attention to very 
red objects which are at the outskirts of the colour distri- 
butions. The red colours of these objects imply they are 
E/SO galaxies at moderate (2 > 1) redshift and consequently 
it is interesting to see if their observed number density is 
consistent with current galaxy formation scenarios. As pre- 



vious authors (iBargcr ct al. 1998 



Moustakas et al. 1997 



Zepf 1997 ; [Cowie et al. 1994 ) who focussed their attentions 



upon these galaxies discussed them in terms of K— selected 
(/ — K) distributions, here we introduce our observations 
in these bandpasses. Following these workers we define an 
"extremely red galaxy" as an object with (/ — K) > 4; from 
Figure ^ we see that at (/ — K) ~ 4, z ~ 1. In this sec- 
tion we will make comparisons with models generated using 
the same set of parameters as the K— selected {B — K) 
distributions discussed previously. In interpreting (/ — K) 
colours one faces some significant differences to [B — K) 
observations; as is apparent from Figure ^ (/ — K) colour 
is relatively insensitive to the amount of star-formation; at 
z ~ 1 the difference in (/ — K) between elliptical galaxies 
with a Scalo IMF and r = 0.5, 1, 2.5 Gyr is 0.5; by compar- 
ison, at this redshift {B — K) ranges from ~ 3 to ~ 5. This 
difference is reflected in the tightness of the K— selected 
(/ — K) distributions plotted in Figure ^. This insensitivity 
makes (/ — K) distributions relatively poor probes of the 
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elliptical star-formation history. However, it is precisely this 
behaviour which makes (/ — K) colours useful in detecting 
high-redshift, evolved ellipticals; Figure ^ shows that until 
a ~ 2 {I — K) colour is approximately proportional to z, and 
is relatively insensitive to variations in r or the IMF. 



Turning to work from the literature, Cowie et al. (1994), 
in a 5.9 arcmin^ survey, detected 13 objects with (I — K) > 4 
to K < 20.9™, and concluded that these galaxies do not 
dominate the faint K > 20"* population. They also con- 
cluded that less that 10% of present day ellipt icals could 
have formed in single-burst events. In contrast, Moustak; 
et al. (l|997|) found 8 galaxies with (/-X) > 4 and < 22 



over a small, 2 arcmin^ area. They also isolated a population 
of objects with blue optical {{V — I) < 2.5) and red near- 
infrared colours and argued that these col ours could not be 
reproduced with Bruzual & Chariot (1993)-type models. Al- 



though we do not have V— band photometry the top mid- 
dle and top left of our Figure ^ roughly corresponds to the 
regions highlighted in the colour-colour plot shown in Fig- 



ure 9 of 



Moustakas et al 



Moustakas et al. 



These "red outlier" objects, as 
describe them, viewed in the context of the 
models discussed here are most lik ely e lliptical galaxies a t 
1 < 2 < 2. Finally, both |Zepf (1997| ) and jBarger et al. (199^ ) 
investig ated extrem ely red objects in the Hubble Deep Field 
(HDF). Zepf (1997), using publicly-available Kitt Peak near- 
infrared imaging data covering the 5 arcmin^ of the HDF 
found that to a 50% completeness lim it of < 22™ ther e 
were ~ 2 objects with (Veoe — K) > 7. Barger et al. (1998), 
in a much larger ~ 60 arcmin^ infrared survey comprising 
both the HDF and the HDF flanking fields found 12 galaxies 
with {I - K) > A a.nd 19 < K < 20. 

Table ^ presents a compilation of observations from 
these papers as well as our current work. We also show the 
predictions of the a; = 3 and Scalo r = 0.5 model (our model 
which most closely approximates a single burst), normalised 
to the total numbers of objects each magnitude slice. Gener- 
ally, in the brighter bin (19 < A" < 20 ) the only sub stantial 
comparison that we can make is with Barger et al. We see 
a factor ^ 3 more objects than they do, although this repre- 
sents only a 2a discrepancy. In the fainter bin 20 < J(" < 22 



our numbers of red objects agree well with [Moustakas et al. 
Although our numbers are higher than Cowie et al. (199^ ) in 



this range there is still no significant discrepancy due prin- 
cipally to our small area. 



Barger et al. compared their observations to the pre- 
dictions of a passively-evolving r = 0.1 Gyr model elliptical 
population, and concluded that the numbers of (7 — K) > 
4 objects observed disagreed with the predictions of this 
model. Our r = 0.5 Gyr distribution (solid dot-dashed 
line) in Figure ^ is quite similar to this model, and in the 
19 < < 20 magnitude slice we can see from Table ^ 
that it does indeed produce a large number of objects with 
{I—K) > 4 whic h are not seen in our observations. We there- 
fore conclude, as Barger et al, did, that the low numbers of 
objects with (I — K) > 4 observed at 19 < A" < 20 disfavour 
non-evolving and passively evolving models (the constraints 
on the models which can be obtained from the distributions 
at fainter magnitudes are less significant). However, we note 
that the x = 3, t = 2.5 Gyr and r — 1.0 Gyr Scalo mod- 
els do correctly reproduce the observed colour distributions. 
This is because both these models contain enough on-going 



star-formation to move galaxy colours sufficiently blueward 
to match the observations. 

Lastly, could a merging model, like those adopted by 



KaufFmann & Chariot (1998), produce these results? It is 



not possible to rule this model out from the redshift data 
and it may even be said that the difficulty that models with 
standard Scalo IMF's have in fitting the K < 19 redshift dis- 
tribution is evidence in favour of the merging scenario. We 
therefore consider whether the K— selected {B — K) and 
(J — K) colour distributions presented here can discriminate 
between a PLE {x — 3) model and a merging model. The 
K— selected (B — K) distributions in Figure ^ can be under- 
stood in terms of either scenario; the red galaxies that are 
missing at 17 < < 22 may either be due to them becom- 
ing blue at high redshift due to continuing star-formation or 
passive evolution - or because they are fainter than expected 
due to de-merging (although our lack of red galaxies even at 
the faintest limits implies either several magnitudes of fad- 
ing or that the pre-merger components are blue). However, 
as we have already noted the (/ — K)-z relation for E/SO 
galaxies is insensitive to changes in r or the IMF. Further- 
more, in the range 18 < ii' < 20 some galaxies with the 
I — K > 4: colour expected of 2 = 1 elliptical galaxies are 
detected (Figure |^. Indeed the numbers of these galaxies 
we see is within a factor of ~ 2 of what is predicted on the 
basis of the Scalo, r = 0.5 Gyr model and in even better 
agreement with the a; = 3, r = 2.5 Gyr models or even a 
Scalo r = 1.0 Gyr model. If the merging model is the ex- 
planation of the large deficiency of red galaxies in {B — K) 
compared to passively evolving models, then it might be ex- 
pected that a similar deficiency should be seen in (/ — K). 
Since the deficiency in (/ — K) is less, then this might be 
taken to be an argument against merging and for the x = 3 
PLE model. 



5 CONCLUSIONS AND SUMMARY 

In this paper we have presented the results of two near- 
infrared surveys to if ~ 20 and if ~ 23 which cover our 
ultra-deep (B ~ 28) optical fields. We draw the following 
conclusions from this work: 

Our K number counts are consistent with the predictions 
of non-evolving models with < go < 0.5 



As previously noted by jSowie, Songaila, fc Hu 



(199e ) jMetcalfe et al. (1996| ), and [Kauffmann fc Chariot 
Tggg) the 1% < K < 19 n{z) of powie, Songaila, fc Hu 



(1996) is also well fitted by non-evolving models. However, 



passively evolving models with a Salpeter/Scalo IMF pre- 
dict too many galaxies with 2 > 1. Dynamical merging is one 
possible solution to reduce the numbers of these galaxies but 
we have shown that a dwarf-dominated IMF for early-types 
could offer an alternative explanation. 

Our K— selected [B — K) colour distributions display a 
strong bluewards trend for galaxies fainter than K ~ 20, 
confirm ing results previously observed in the shallower sur- 
veys of Gardner, Cowie, fc Wainscoat (1993 ). 

At brighter magnitudes (K < 20 mag) our K— selected 
{B — K) distributions indicate a deficiency of red, early-type 
galaxies at 2 ~ 1 compared to the predictions of passively 
evolving models. This implies either a PLE model where 
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Table 5. Numbers of objects with (/ — K) > 4 per arcmin-^ with itltr errors. The survey area in arcmin^ is shown in parentheses. In the 
deep survey (20 < K < 22) there are 7 objects which are undetected in / and therefore only lower limits can be placed on their colours. 
They are included in their respective bins. The number quoted from Moustakas et al. (1997) covers the magnitude range 20 < K < 22. 



Author 

Model (x = 3) 

Model (scalo, t = 1.0 Gyr) 

Model (scalo, r = 0.5 Gyr) 

This work 



Barger et al. (1998) 



Cowic ct al. (1994) 



Moustakas ct al. (1997) 



19 < K 


< 20 


20 < i^: < 21 


0.3 




1.4 


0.6 




1.4 


1.1 




3.1 


5+"-^ 
'^■°-o.i 


(47) 


3.9tll (1.8) 


2+°-^ 


(62) 


1.2t°i (7.8) 


3+°" 


(7.1) 


(5-9) 






4.0+' 



21 < if < 22 
1.8 
1.2 
2.8 

1.1 ' ' 



hi. 5 
'0.7 



0.3: 

(2) 



(1.8) 
(5.9) 



star-formation continues at a low level after an initial burst 
or dynamical merging. 

At fainter magnitudes (20 < K < 22) the continuing blue- 
wards trend observed in {B — K) can be explained purely 
in terms of passively evolving PLE models with no need to 
invoke any additional mechanisms. 

Our observed numbers of (/ — K) > 4 galaxies &t K ^ 2Q 
are lower than the predictions of passively evolving models 
or PLE models with a low level of continuing star-formation, 
suggesting that at least part of the larger deficiency observed 
in [B — K) at 7^ ~ 20 may be due to star-formation rather 
than dynamical merging. 

In the range 19 < AT < 20, where our statistical uncer- 
tainties are lowest, we detect 0.5±0.1 red galaxies arcmin^^ 
with (J — K) > 4. We see a factor ~ 3 more objects than 



Barger et al. do although this represents only a 2a discrep- 
ancy. The PLE models discussed here suggest these galaxies 
will have redshifts 1 < z < 2. The numbers of these galax- 
ies are consistent with the predictions of PLE models with 
small amounts of ongoing star-formation. 



6 FUTURE WORK 

We have recently completed a H~ band survey of the Her- 
schel field using the 1024 x 1024 Rockwell HgCdTe array at 
the Calar Alto 3.5M telescope. This survey is complete to 
H < 22.6'" and covers the entire area of the Herschel Deep 
Field. Analysis of this dataset will be presented in future 
papers. 
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